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An elementary pedagogical derivation of the Lense-Thirring precession 
is given based on the use of Hamilton vector. The Hamilton vector is an 
extra constant of motion of the Kepler/ Coulomb problem related simply to 
the more popular Runge-Lenz vector. When a velocity-dependent Lorcntz- 
like gravitomagnetic force is present, the Hamilton vector, as well as the 
canonical orbital momentum are no longer conserved and begin to precess. 
It is easy to calculate their precession rates, which are related to the Lense- 
Thirring precession of the orbit. 

PACS numbers: 45.20.D-, 91.10.Sp, 04.25.Nx 

1. Introduction 

Every relativistic theory of gravitation must include a Lorentz-like force 
induced by a magnetic- type component of the gravitational field; the general 
theory of relativity by Einstein does so. More precisely, in its weak-field and 
slow-motion approximation the highly nonlinear Einstein field equations get 
linearized, thus resembling the linear equations of the Maxwellian electro- 
magnetism. As a consequence, a magnetic-type component of the gravita- 
tional field appears, induced by the off-diagonal components goi, i = 1, 2, 3 of 
the spacetime metric tensor [1] . The role of the electric currents is played by 
mass-energy currents: in the case of a slowly rotating mass, a test particle 
geodesically moving far from it is acted upon by a non-central velocity- 
dependent Lorentz-like force. For other effects induced by the gravitoma- 
gentic field on the motion of gyroscopes, test particles, moving clocks and 
atoms, light rays see, for example, [2]. 

Recent years have seen increasing efforts, both from theoretical and ob- 
servational points of view, towards a better comprehension of the post- 
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Newtonian gravitomagnetic field: for a comprehensive recent overview see, 
for example, [3]. Thus, we feel it is not worthless to offer to the reader an 
elementary derivation of the precessional effects induced by the gravitomag- 
netic field on some Keplerian orbital elements of a test-particle, i.e. the so 
called Lense-Thirring effect [1] , although a recent historical analysis [5] sug- 
gests that it should be more appropriately named Einstein-Thirring-Lense 
effect. For other derivations of such an effect see, for example, [6] [71 IHl l9l [TO] . 

Our presentation is aimed for the first year physics students with limited 
experience in mathematics and theoretical methods of physics. We, there- 
fore, give somewhat detailed exposition which an experienced physicist may 
find unnecessary long, but as our teaching experience shows such an ex- 
position is helpful and necessary for newcomers in the field. We hope the 
material presented here will help students which just begin their physics ed- 
ucation not only master some simple mathematical methods learned during 
the first year mechanics course but also feel the beauty of advanced topics 
which they will learn in detail later on. 

2. Gravitomagnetism 

Nowadays it is widely accepted that the correct theory of gravitation is 
provided by the Einstein's general theory of relativity [II] (for an excellent 
introduction for beginners see |12j : a classic textbook is p3]). However, the 
general theory of relativity is a highly non-linear theory. The spacetime in 
it is not merely a static background for physical processes; it is dynamic 
and affected by any contribution to the energy-momentum tensor of the 
system under investigation. Such a tensor enters in a prescribed manner 
the Einstein's equations which determine the ten components of the metric 
tensor. This metric tensor enters by itself the equations of motion of the 
system under study. There is no hope for a neophyte to struggle his /her way 
into this impossibly tangled up mess of non-linear jungles "unless months 
of study on the specialized terminology, procedures, and conventions of the 
general relativity theorist have been completed" |14j . 

However, virtually the large part of all the performed and/or proposed 
tests, aiming to observationally scrutinize the post-Newtonian regime, deal 
with weak fields and non-relativistic velocities. Therefore, we may expect 
that the full machinery of general relativity is not necessarily required to 
estimate the relevant post-Newtonian effects. Some linearized version of 
the theory will suffice to do the job much faster and physically in a more 
transparent manner [SI [T^l [TdI [IB] . 

Although, in general, spacetime is dynamical and there is no natural 
way to split it into space plus time, for stationary spacetimes, like the one 
around the Earth, stationarity dictates a preferred way of how this splitting 
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can be performed. There is no approximation in such slicing of spacetime 
into three-dimensional space plus one-dimensional time. It is just a new 
mathematical language convenient in stationary situations |16| . 

Under 3 + 1 slicing, the spacetime metric tensor g„ v naturally decom- 
poses into several parts. In the case of weak gravity and non-relativistic 
velocities, this decomposition allows one to set up a remarkable analogy 
with electromagnetism [21 [HI HU [16] . In particular, gravitational analogies 
of the electromagnetic scalar and vector potentials, <& and A respectively, are 
determined by the time-time and time-space components of the spacetime 
metric: ^ 

where c is the speed of light. Gravitoelectric field E and gravitomagnetic 
field H are related to these potentials (in the Lorentz gauge [2]) in the 
usual way (up to an extra factor 4 which also appears in the Lorentz gauge 
condition): 

1 dA 

E = -V$-— — , H = VxA 

4c dt 

and satisfy a gravitational analog of Maxwell's equations [2[ [TH[ [T7] 
V • E = -AttG P , V • H = 
Vxl = 0, V x H = 4 



4ttG P ^ + 1 ® E 
c c dt 



(1) 



Similarity with Maxwell's electromagnetic equations is apparent. The role 
of charge density is played by mass density, p, times Newton's gravitation 
constant, G. The mass current density, Gpv, with v as the velocity of the 
source mass, plays the role of the charge current density. 

However, there are several important differences [16 J. Gravity is medi- 
ated by a spin- two field and is attractive. In contrast, electromagnetism 
is mediated by a spin-one field and can be both attractive and repulsive. 
This difference leads to the extra minus signs in the source terms in ([1]). 
Another remnant of the tensor character of gravity is an extra factor, 4, in 
the equation for V x H. No gravitational analog of the Faraday induction 
is present in (pQ). However, this is an artifact of restricting to the only first 
order terms in v/c of the gravitating masses in the above equations. At the 
second order, {—l/c){dH/dt) term reappears in the equation for V x E, 
but also some non-Maxwell-like terms are introduced elsewhere in the field 
equations [151 02] • 

Once the gravitoelectric and gravitomagnetic fields are known, the force 
acting on a small test body of mass m is given by the formula which is 
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analogous to the Lorentz force law 

F = mE H v x H. 

c 

More precisely, there is still another difference from electromagnetism as 
far as the equation of motion is concerned. So far nothing was said about 
the space-space components, gij, of the metric tensor. In general, they 
do not correspond to Euclidean space but to that of curved space. Just 
the space curvature effects are responsible for the classic general relativistic 
predictions for Mercury's perihelion precession and bending of light rays. In 
the presence of space curvature, we have the curvilinear equation of motion 

ld(Tvi) l(dg jk \ jk „ 1_ ldHi 

i » v> v * = Ei + - (v x H)i — , (2) 



where 



T dt 2\dx i J c v ' c dt 

2* v'Y 1 ' 2 



1 



c 2 c 2 



and, as usual, summation over repeated indexes is assumed. 

In this paper, we are interested in tiny secular effects of Earth's rotation 
on a satellite orbit. Therefore, we will neglect the space curvature effects, as 
well as terms of the order of (v/c) 2 other than gravitomagnetic. Although 
these curvature effects are even more prominent than the effects we are 
interested in, they are well known. Therefore, for our goals, equation ^) 
takes the form 

dv a _ m 

— = — j n+ ~ 

dt r z c 

where r = rn is the radius-vector of the satellite of mass m, a = GmM, 
where M is the Earth's mass, and H is the gravitomagnetic field caused by 
Earth's rotation. 

The gravito-electromagnetic analogy implied by (pQ) allows us to find H 
|17| . It is well known that the magnetic moment 



ix _, ill/ 

m— = 2 nJ[ vxH, (6) 



related to the electric current density j, creates a dipole magnetic field 

3n (n • p) — fl 



H 



Equations ([T]) indicate that the gravitational analog of the electric current 
density is —4Gpv. Therefore the gravitational analog of the magnetic mo- 
ment is 

fl = -AG — [ p[rxv]dV = -2G -, 
2c J c 
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where S = f[r x v] pdV is the rotating body's proper angular momentum. 
Thus the gravitomagnetic field of the Earth is given by 

n 2G S-3n(n-S) 

H = — -s • ( 4 ) 



3. Larmor precession 

It is instructive first to consider in detail a simpler case of a non- 
relativistic Coulomb atom slightly perturbed by a weak and uniform con- 
stant magnetic field H. The equation of motion is 

l, ol _ e _, - 
mv = — ^ n H — v x H, (5) 
r 2 c 

where now m and e are electron mass and charge respectively, and a is the 
strength of the corresponding electromagnetic coupling between the electron 
and the atom nucleus. 

In the absence of perturbation, H = 0, equation (|5|) describes the elec- 
tron motion on the ellipse. It is convenient to characterize the orientation 
of this ellipse by the angular momentum vector L = fx p, which is perpen- 
dicular to the orbit plane, and the Runge-Lenz vector 

A = v x L — an, (6) 

which is directed towards perihelion of the orbit and equals in magnitude 
to the orbit eccentricity times electromagnetic coupling a. The Runge- 
Lenz vector is an extra constant of motion of the Coulomb problem and its 
existence is related to the so called hidden symmetry of the problem [18\ [T9] 
which makes the Kepler/Coulomb problem an interesting testing ground 
for various algebraic and geometric methods [201 ED Note that first 
integrals for the Kepler Problem, equivalents of the Runge-Lenz vector, were 
first obtained by Ermanno and Bernoulli many years before Gibbs devised 
his vectorial notation [23, 24, 25J. "The lack of recognition of Ermanno and 
Bernoulli parallels that of Aristarchos and makes one wonder if there is some 
sinister influence for those who are initiators in the field of the motions of 
the planets" [23] . 

The Runge-Lenz vector A and the angular momentum vector L are 
mutually perpendicular. Therefore we can find the third vector u such that 

A = u x L. (7) 

This vector is called the Hamilton vector [261 EH ESI EH] and it has the form 



u = v - — e v , 



(8) 



G 
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where 

= — - L x n = I x ft, with I = — , (9) 

is the unit vector in the direction of the polar angle ip in the orbit plane. 
Sometimes it is more convenient to characterize the orbit orientation by 
vectors u and L. 

Now let us return to © with the non-zero magnetic field H. As the 
Lorentz force, in general, has a component outside the unperturbed orbital 
plane, it is evident that the motion will no longer be planar. However, 
for weak magnetic fields, the Lorentz force component in the orbital plane 
is much smaller than the binding Coulomb force. Therefore, it can cause 
only small perturbations of the elliptical orbit. The Lorentz force compo- 
nent perpendicular to the orbital plane is expected to cause this plane to 
turn slowly around the Coulomb field center. This intuitive picture of the 
perturbed motion suggests us to decompose the electron velocity in such a 
way 

v = v' + ttxr. (10) 

Here v' is the electron velocity relative to the instantaneous orbital plane, 
and the second term, Vt x r, is due to slow revolution of this plane with a 
small angular velocity ft. 

If we want the angular momentum L to be perpendicular to the instant 
orbital plane, we should replace v by v' in its definition: 

L = rxP, P = mv'. (11) 

Analogously, the appropriately generalized Hamilton vector is 

u = v — — e<p = T e ¥" ( 12 ) 

L m L 

where e v is still given by equation ([9|). 

The Runge-Lenz vector, defined as earlier by ([7]), is no longer conserved. 
However, at instances when the particle is at perihelion of the instant ellipse, 
the Runge-Lentz vector still points towards the perihelion, as the following 
simple argument shows [30]. At the perihelion of the instant ellipse, the in- 
trinsic velocity v' has no radial components and, therefore, is perpendicular 
to the radius vector f. Then (|12p shows that the Hamilton vector is also 
perpendicular to r and, according to ([7]), the Runge-Lenz vector A will be 
parallel to r, that is pointing towards the perihelion. 

If our intuitive picture of the orbit plane precession is correct, we should 
have 

L = nxL. (13) 
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On the other hand, 

L = fxP + fxP. 
As for the velocity v = f , we have from (|10p 

P 

r = h Q X r . 

m 

While, up to the first order terms in H and 0, 



d i 

P = — [mv — m\l x r 
dt 



a e -> -> -> ^ 

■—^n-\ P x H — 0, x P — ml] x f. 

r z mc 



Substituting (TT5]) and (TT6D in CED> we g et 

L = (Qxr)xP + — r x (P x H) — r x ($1 x P) — mf x (!1 x r) 

mc 

or, after using (0 x r) x P = —r X (0 X P) — X (P X r), 



L = l]xL-fx 



2f2 H P) x P 

mc 



mf x(Oxr). 



As we see, to get (fT3j) . it is sufficient to take 





6 P. 



2mc 



Now let us calculate (note that r = v ■ n) 



Substituting 
and using 

we get 
Analogously, 



v f(v ■ n) 1 
n = ^ — = — [n x v) x n. 



P 



v = h x r 

m 



(n x P) x n = — eL, [n x (0 x n)l x n = (l x n, 
r 



n = Q x n H ^ eL. 



L 



e (/3 = /xn + /xn = (J7x/)xn + /x(r2xn)H ^ ' x e 



mr 
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But (Q x T) x n + I x (O x n) = — f2 x (n x I) = x and I x = — n. 
Therefore, 

ei» = x 4 ^ ft. (19) 

At last, by using ([5j) for u, we get 

V « ^ e ^ ^ _ a ^ 

f = n H uxii-lixu = Ux« n. 

mr 2 mc mr 2 

But, up to first order terms in f2, f2 x v m $7 x v'. Therefore, 

v &ttxv' ^- n. (20) 

mr 

Having ([20]) and (fT9|) at hand, it is easy to find 

-> a -* -* 

u^Qxv' — — £lxe ip = Qxu. 

Therefore, the Hamilton vector precesses with the same angular velocity 
(|17p as the angular momentum, and this angular velocity can be considered 
as related to the precession of the instantaneous ellipse as the whole. 

As the final remark, note that A = x r is the vector-potential cor- 
responding to the uniform magnetic field H . Therefore 

-* -* 6 -* 

P = mv — m!l x f = mv H — A 

c 

is just the canonical momentum of the non-relativistic electron in the mag- 
netic field H. 

4. Lense-Thirring precession 

It is well known [31] that the vector-potential created by a magnetic 
moment jl is 

*W = ^ (2D 

Therefore, by making the change jl — > —^-L' in (|2ip . we get the gravita- 
tional analog of the vector potential for the gravitomagnetic field ([!]) of the 
Earth with angular momentum L'\ 

1C 

A = -^rxL'. (22) 



cr 3 



To recast the canonical momentum, 



m 

P = mv -\ A, 

c 
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in the form P = mv — mil, x f, which corresponds to the decomposition of 
velocity (fT0j) . we could take 

^ = % L '- ( 23 ) 
Then the gravitomagnetic field ([4]) can be rewritten as 

— = n-3n(n-n). (24) 
c 

By using (j24"|) and equation of motion ([3]) , we get 

-> ot -> -> -* 

P = — g ^ + 2m C x — 3m (0 • n) u x n — m x r. (25) 

From ([23]) . we find 

Q. = —3— (v ■ ft), 
r 

and (1251) takes the form 



P = — g ™ + 2m u x — 3m [u (f2 ■ n) — f2 (v ■ ft)) x n. 

However, 

[v (Q • ft) — Q (v ■ fi)] x n = [ft x (u x fi)] xn = t/xf2 — ft ft - (v X Q), 
and we end up with the equation 

P = — -n + m!]xi) + 3mnii'(iJx(]). (26) 



Up to the first order in the small parameter £1, (|26p can be rewritten as 

Pra— ^n + QxP + 3iifi'(nxP) (27) 
Now it is easy to find 

L = ( — - - A ) xP + rx(OxP) = -(f x Q) x P + r x (O x P). (28) 
\m c / 

But r x (tt x P) + P x (f x Q) = -O x (P x r) and ([25]) indicates that the 
canonical angular momentum vector L really precesses with angular velocity 
O: 

L = nx L. (29) 



10 



LT' Hamilton printed on May 26, 2009 



We are interested in the secular changes of the orbital parameters. Therefore 
it makes sense to average Q in (|29p over fast orbital motion: 



2C 1 

n^ < n > =—L'<- E >. (30) 

For the desired accuracy, we can average 1/r 3 in dUD over the unperturbed 
orbit 

- = 1 + e cos (p, 
r 

where p and e are the semi-latus rectum and eccentricity of the orbit. But 
for the unperturbed orbit 

2 

dt = —— dip, (31) 
1j 



T 2tt 

1 1 f dt m f, . 2irm 



and we get 

f dt m /\ 

" r 3 T J = ~LTp J C ° S ~ LTp ' 

o o 

Integrating ([BT]) over the complete orbital period T, we get 

LT 



= 2S = 2vraVl - e 2 , 

777 

where S 1 is the area of the ellipse and a is its semi-major axis. On the other 
hand, p = a(l — e 2 ) and we finally get 

1 _ 1 

< ^3 >- a 3 (1 _ e 2)3/2- ( 32 ) 

Therefore, the averaged angular velocity of the precession is 

-> 2G 
<fi>= c 2 a3(l- e 2 )3/2 L '- ^ 



It is not difficult to check that relations (I18j) and (j 19[) remain valid if 
Q, in these formulas is given by ([23]) . Then, after using ([27|) and (fl9|) . we 
easily find 

77 2G 77 

u^Clxu + S Q - L = tixu + 3 — k—tL'-L. (34) 

mr mc z r 4 

Therefore, the Hamilton vector does not precess with angular velocity fl. 
However, we should average ()34|) over the fast orbital motion. The first term 
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simply gives <f2xu> = <f2> x u because u is a slowly changing vector 
and we can assume that it does not change over time scales comparable to 
the orbital period T. As for the second term, we can use ft = cos ip i + sin p j 
and get 

2tt 

n m f. . , 2 , 27rme . . 

< r 3 >= TLp 1 J ( COS(fil + sm( PJ)\ l + ecos^j d( P = YZ^ 1 ' ^ ' 

However, 



o 



-* -f j* u L e 1 m 

1 — j X /i — — X — . — — — — ^ 

it L pnL pa L 2 

and (1351) takes the form 



< -r >=< ^ > -ru x L. (36) 



Therefore, after averaging, (JMJ) changes to 

u Q LT x u, (37) 

where 

of 

aLT= c^(l-e^ '- 3l{l - Z ' )] - (38) 

As we see, the secular precession of the elliptic orbit contains two terms: the 
precession of the orbital plane with the angular velocity < 0, > around the 
central body's angular momentum L' , and the precession within the orbital 
plane with the angular velocity £Ilt — < ^ > around the angular momen- 
tum L. The magnitudes of the angular momentum L and the Hamilton 
vector u, and therefore, the orbital parameters such as the eccentricity and 
the semi-major axis, remain unchanged to the first order of the perturbation 
theory. 



5. Concluding remarks, Skovoroda's principle and all that 

Lev Borisovich Okun cites [32j eighteenth century Ukrainian philosopher 
Grigory Skovoroda as the author of the remarkable principle, which Okun's 
PhD adviser Isaak Yakovlevich Pomeranchuk used to quote: "Thanks God: 
All what is relevant is simple, all what is not simple is not relevant." 
The formal analogy between weak field low velocity general relativity and 
Maxwellian electrodynamics is a simple and elegant way to illuminate a 
whole class of interesting physical phenomena dubbed gravitomagnetism. 
Lense-Thirring precession is one such example. However, it should be kept 
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in mind that the analogy is only formal and sometimes can lead to strange 
and erroneous conclusions if we forgot about rather strong limitations under 
which the approximation underlying the analogy is valid [33J . 

Frame dragging is maybe more appropriate interpretation of the Lense- 
Thirring effect not restricted to the weak field limit. However, "whereas 
'frame dragging' is a very catchy appellation" [33] its meaning is not easy 
to explain to introductory level students. Therefore, for them and not only 
explaining the Lense-Thirring effect by analogy with Larmor precession re- 
mains a relevant and simple option. 

The use of Hamilton or Runge-Lenz vector greatly simplifies the dis- 
cussion. However, a great deal of vector algebra is still needed in either 
cases, as the previous chapters illustrate. This is in contrast with situation 
in perihelion precession under central force perturbations where the use of 
the Hamilton vector trivializes the problem [35]. In fact, the amount of 
vector algebra is a price we should pay for our desire to keep the approach 
elementary and accessible in introductory mechanics course. If more ad- 
vanced background in analytical mechanics is assumed, much more elegant 
and technically simple way is provided by the use of the Poisson brackets 
[36] or, alternatively, by Hamilton equations of motion [37] . 

In terms of the standard Poisson brackets, we have 



The Poisson brackets are easy to calculate by using the Leibniz rule 
{fi9h} = {f,g}h + g{f,h} and the fundamental Poisson brackets 



where f(r,p) is any scalar function of its arguments. It follows from (|40p 
that 




{n,rj} = 0, 

{Li, Tj} = Eijk Vk 




{Li,f(r,p)} = 0, 



(40) 




(41) 



Therefore 



which is equivalent to (|29p . 
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Besides, for any vector of the form B = fi(r,p)f+ f2(r,p) P, with any 
scalar functions /i and and in particular for vectors e v and u, we will 
have 

{Li, Bj} = e^k Bk, 

and we get immediately 

1 3 n _ _ 
■hi = {ui, Lj}Vtj H {Pi, Qj}Lj = e ijk u k Qj H Q ■ L, 

which is equivalent to (|34j) ■ 

As we see, in combination with some background in Hamiltonian me- 
chanics, the use of Hamilton vector again makes the exposition rather trivial, 
in complete agreement with the Skovoroda's principle. 

In our discussions canonical variables played an important role, as they 
are in accord with the intuitive physical picture of precessing orbital plane. 
However, the use of canonical variables are not mandatory for the pertur- 
bation theory discussion of the Lense-Thirring effect. The standard per- 
turbation theory in celestial mechanics is based on variation-of-parameters 
method which has some inherent gauge freedom [38J. Indeed, the unper- 
turbed Keplerian ellipse is determined by six parameters 

r = f(C 1 ,C 2 ,C 3 ,C 4 ,C 5 ,C 6 ,t). (42) 

In the role of these parameters one can take, for example, three Euler an- 
gles which determine the orbit plane orientation, the semi-major axis and 
eccentricity of the orbit, and the so called mean anomaly at an epoch, which 
determines the initial position of the body. Another obvious choice is the 
initial values of the position and velocity vectors of the body, and many 
other sets (Delaunay, Poincare, Jacobi, Hill) can be found in the literature 

M- 

To solve the perturbed equation 

•• ct — 
mr = -n + F, (43) 

with F as a small perturbation, by the variation-of-parameters method, 
one can assume that the solution still has the form (|42p but Cj are no 
longer constant. However, three scalar equations (|43[) are not sufficient to 
determine six unknown functions Cj(i). We need three auxiliary conditions 
on them and the freedom in choosing of these auxiliary conditions is just 
the gauge freedom mentioned above. Usually the gauge fixing is achieved 
by the Lagrange constraint 
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Under this condition 




That is, the velocity is the same function of the parameters Cj as in the 
absence of perturbations and the instantaneous ellipse, determined by these 
parameters, is tangent to the real trajectory. Correspondingly, the instan- 
taneous orbital parameters Cj so determined are called osculating. 

However, not all choices of orbital parameters Cj are compatible to the 
Lagrange constraint (|44|) . For example, in previous chapters we have choose 
the canonical angular momentum and the Hamilton vector to characterize 
the orbit shape and orientation. When velocity dependent perturbations are 
present, such parameters are not osculating, as is evident from the ascribed 
physical picture of the combined motion. 

The canonical perturbation theory in celestial mechanics was developed 
in [8]. Again the corresponding orbital elements are not osculating if the 
perturbation depends on velocity. In fact this is a general property: under 
velocity-dependent disturbances, canonicity and osculation are not compat- 
ible [31]. 

The problem of orbit perturbations can be solved in any gauge. At 
that some gauges are more convenient, because if the gauge corresponds to 
the real physical picture of the perturbed motion, the formalism simplifies. 
Some examples are given in |8j when the standard perturbation theory with 
osculating elements is cumbersome while the canonical perturbation theory 
is more elegant. 

The Lense-Thirring effect can be considered either in the canonical per- 
turbation theory, or with the osculating orbital elements [8] , or in any other 
convenient gauge. At that the corresponding orbital elements can differ con- 
siderably. Of course, the real physical quantities, like position and velocity, 
do not depend on the gauge used after the initial values are accounted cor- 
rectly [8]. However, sometimes the orbital elements, like inclination angle, 
are also considered as physically real. In such cases care should be taken to 
relate the measured quantities to the orbital parameters used in the pertur- 
bation theory. Without this care a confusion can arise when two different 
mathematically correct approaches give seemingly different results which 
are in fact equivalent [6]. 
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